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ABSTRACT 

We investigate how the presence of a non-thermal tail beyond a Maxwellian electron 
distribution affects the synchrotron process as well as Comptonization in plasmas with 
parameters typical for accretion flows onto black holes. We find that the presence of 
the tail can significantly increase the net (after accounting for self-absorption) cyclo- 
synchrotron emission of the plasma, which emission then provides seed photons for 
Compton upscattering. Thus, the luminosity in the thermally-Comptonized spectrum 
is enhanced as well. The importance of these effects increases with both increasing 
Eddington ratio and the black hole mass. The enhancement of the Comptonized syn- 
chrotron luminosity can be as large as by factors of ~ 10^ and ~ 10^ for stellar and 
supermassive black holes, respectively, when the energy content in the non-thermal 
tail is 1 per cent. 

The presence of the tail only weakly hardens the thermal Comptonization spec- 
trum but it leads to formation of a high-energy tail beyond the thermal cut-off, which 
two effects are independent of the nature of the seed photons. Since observations of 
high-energy tails in Comptonization spectra can constrain the non-thermal tails in the 
electron distribution and thus the Comptonized synchrotron luminosity, they provide 
upper limits on the strength of magnetic fields in accretion fiows. In particular, the 
measurement of an MeV tail in the hard state of Cyg X-1 by McConnell et al. implies 
the magnetic field strength in this source to be at most an order of magnitude below 
equipartition. 

Key words: accretion, accretion discs - gamma-rays: theory - radiation mechanisms: 
thermal - X-rays: galaxies ~ X-rays: stars. 



1 INTRODUCTION 

Power-law-like X-ray spectra extending up to soft 7-rays 
(hereafter X7) are a common feature of the emission from 
accreting black holes in active galactic nuclei (AGNs) and 
black-hole binaries (BHBs). It is generally accepted that 
those spectra are produced by the Comptonization process 
taking place in optically thin accreting plasmas with typical 
electron energies of ~ 100 keV. Similar spectra are also ob- 
served from weakly-magnetized accreting neutron stars (e.g. 
Barret et al. 2000), where they can originate either in the 
accretion flow or in a boundary layer near the neutron star 
surface. 

The X7 spectra of Seyfert Is and X-ray binaries in their 
hard (low) states can be modelled by Comptonization on 
thermal electrons, see e.g. Zdziarski (1999). However, both 
observations and theoretical considerations suggest that a 
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small non-thermal component may be present in the electron 
distribution, giving rise to a weak non-thermal high-energy 
tail in the observed Comptonization spectrum. Such a tail 
above ~ 1 MeV has been detected in the spectrum of Cyg 
X-1 in the hard state by the COMPTEL detector aboard 
CGRO (McConnell et al. 1994, 2000a). In other BHBs and 
in Seyferts, the existing upper limits are compatible with 
the presence of weak non-thermal tails (e.g. Zdziarski et al. 
1998; Gondek et al. 1996; Johnson et al. 1997). 

A non-thermal tail in the electron distribution in com- 
pact X7 sources can arise due to many difi'erent physical 
mechanisms. One such mechanism is acceleration in the pro- 
cess of dissipation of magnetic field (similarly to the case of 
the solar corona) in an optically-thin accretion flow or in 
active coronal regions above an optically thick disc. Then, 
the energy loss and subsequent thermalization of the accel- 
erated electrons leads to a steady-state electron distribution 
consisting of a Maxwellian and a non-thermal tail. Apart 
from Coulomb interactions, an effective thermalizing pro- 
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cess is synchrotron self-absorption (Ghisellini, Guilbert & 
Svensson 1988). 

A strong support for the occurrence of non-thermal pro- 
cesses in black-hole accretion flows comes from observations 
of X7 spectra of black-hole binaries in the soft state. In 
that case, Compton scattering by non-thermal electrons ap- 
pears to be the dominant process (Poutanen & Coppi 1998; 
Gierliiiski et al. 1999). In the well-studied case of Cyg X- 
1, the power-law-like spectrum from this process becomes 
dominant at ^ 10 keV (Gierlihski et al. 1999) and it extends 
to at least ~ 10 MeV (McConnell et al. 2000b). 

An important issue concerns the source of soft photons 
undergoing Comptonization in the accreting plasma. A nat- 
ural candidate is thermal emission from an optically thick 
accretion disc, either below the corona or outside the re- 
gion of optically thin accretion flow. Ifowever, as the pres- 
ence of magnetic fields of the order of that corresponding to 
equipartition with gas pressure (i.e. ~ 10^ G in AGNs and 
~ 10^ G in BHBs) is expected, one has to take into account 
Comptonization of synchrotron photons (hereafter the CS 
process) when calculating the emergent spectra. In models 
of advection-dominated accretion flows, the CS process may 
be the main radiative process in the innermost part of the 
accretion flow, where most of the available gravitational en- 
ergy is liberated (e.g. Narayan & Yi 1995). 

Wardzmski & Zdziarski (2000, hereafter WZOO) have 
investigated the efficiency of the thermal CS process under 
physical conditions relevant for accreting black holes and 
found that it can dominate only in low-luminosity AGNs 
and in luminous stellar X-ray sources with the hardest spec- 
tra. In the latter case, however, the presence of a correlation 
between the spectral index and the strength of X-ray reflec- 
tion (Zdziarski, Lubihski & Smith 1999; Gilfanov, Revnivt- 
sev & Churazov 1999; Lubihski & Zdziarski 2001) appears 
to imply that blackbody emission of cold media is the main 
source of seed photons for thermal Comptonization, and the 
CS process is negligible. 

In this work, we extend our previous results to the case 
of a hybrid, quasi-thermal plasma and investigate the influ- 
ence of a weak non-thermal component in the electron dis- 
tribution on the efficiency of the CS process. In Section^, we 
discuss the shape of the electron distribution function result- 
ing from acceleration processes operating in accretion flows. 
The influence of non-thermal tails on the CS emission of 
the plasma is discussed in Section ^ We then investigate in 
detail synchrotron emission and Comptonization spectrum 
of hybrid electrons in Sections ^ and respectively. This 
allows us to estimate the efficiency of the CS process in hy- 
brid plasmas (Section ^.2| ), and those results are applied to 
astrophysical black-hole sources in Section 5.3, Finally, we 
summarize our results in Section ^ 



2 HYBRID ELECTRON DISTRIBUTIONS 

In accretion flows, hybrid electron distributions can be pro- 
duced e.g. by stochastic acceleration of thermal electrons 
by interaction with plasma instabilities, e.g. with whistlers 
via gyroresonance (as applied to the case of Cyg X-1 by Li, 
Kusunose & Liang 1996) or with fast-mode Alfven waves 
via transit-time damping (Li & Miller 1997). Such process 
can result from dissipation of magnetic field in an accretion 



disc corona (where magnetic fields are raised from the disc 
by buoyant forces, Galeev, Rosner & Vaiana 1979; Haardt, 
Maraschi & Ghisellini 1994) or in an optically thin accretion 
flow (Bisnovatyi-Kogan & Lovelace 1997), provided a signif- 
icant fraction of the released energy generates turbulence. 

Li & Miller (1997) considered parameters typical for a 
plasma in an accretion disc corona around a stellar black 
hole. They have shown, by solving the Fokker-Planck equa- 
tion, that electrons can be effectively accelerated from the 
thermal background to relativistic energies, so that in the 
stationary state a non-thermal tail of a power-law-like shape 
develops above certain Lorentz gamma factor 7nth ~ 2 
(which value is determined by the efficiency of thermaliza- 
tion processes, see Section ^,) and the power-law index, p, 
depends strongly on the plasma optical depth. The power- 
law component is cut off at a 75 ~ 10 due to radiative losses 
overcoming acceleration. Dermer, Miller & Li (1996) pre- 
dict somewhat larger cut-off energies in the case of electron 
acceleration by whistlers and in a plasma with parameters 
relevant for AGNs. 

Similar shape of steady-state electron distribution is 
also predicted from thermalization by synchrotron self- 
absorption and Coulomb interaction of an injected power- 
law electron distribution (produced by any acceleration 
mechanism), as shown by Ghisellini et al. (1988). 

The exact shape of the non-thermal tail depends 
strongly on the plasma parameters and the (uncertain) 
mechanism of magnetic field dissipation and electron accel- 
eration. We shall therefore consider a generic electron distri- 
bution function of the form qualitatively reproducing that 
predicted by detailed modelling of electron acceleration, 

"^^'^^ ^ I n-(e,7„«0(^)~'exp(2^), 7 > ,J'^ 

Here O = kT/m^c? is the dimensionless electron tempera- 
ture and n]^{Q,^) is the relativistic Maxwellian normalized 
by the relation, no = no(7)d7, where no is the total 
electron density. This hybrid electron distribution consists 
thus of a thermal part below 7nth and an e-folded power-law 
above it (see also Section ^). 

We note that the cases p 3 and p > 3 are qualitatively 
different. In the former, most of the non-thermal synchrotron 
luminosity is produced by electrons with 7 ^ 7f and the non- 
thermal synchrotron spectrum is hard. In the latter case, 
the luminosity is produced mostly by low-energy electrons, 
7 ~ 7nth, and the spectrum is soft. Observations of Cyg X- 
1 (Gierliriski et al. 1999; McConnell et al. 2000a) suggest 
p > 3. In the following we assume p ^ 3. 

The degree to which a hybrid distribution is non- 
thermal can be measured by the quantity, (7nth — 1)/©) 
which shows how far in the Maxwellian tail the non-thermal 
power law starts. Another such quantity is the ratio of the 
energy densities in the non-thermal and thermal parts of the 
electron distribution, 5, which, for mildly relativistic temper- 
atures and 7nth — 1 3> O can be approximated as 
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Figure 1. The relation between jnth E'nd <5 for different plasma 
parameters. 



where K2 is a modified Bessel function and we employed an 
approximation for the thermal energy density from Gammie 
& Popham (1998). In the above derivation (as well as in all 
analytical approximations below) we neglected, for simplic- 
ity, the exponential cut-off in the distribution (0). Typical 
relations between S and (7nth — I)/© ai'e shown in Fig. |l|. 

Since the exponential term in equation ^ varies fastest, 
we have an approximate dependence, 

7„ft~C-ein5, (3) 

where C ~ 1 is slowly increasing with increasing O and 7nth. 
We then see that (7nth — 1)/© only weakly depends on O 
and the shape of the tail for a given S, see Fig. m 




V [Hz] 



3 GENERAL FORM OF EMISSION FROM A 
HYBRID ELECTRON DISTRIBUTION 

In accreting black holes, synchrotron emission from thermal 
plasmas optically thin to scattering is usually strongly self- 
absorbed up to a turnover frequency, Vt, below which the 
plasma is optically thick to absorption, 

Q,,i?=l, (4) 

where is the absorption coefficient and R is the charac- 
teristic size of the source. Hereafter, we assume that Vt is 
determined only by synchrotron absorption, i.e. other radia- 
tive processes are negligible at v Ut. Typically, accretion 
flow models predict (assuming equipartition magnetic fields) 
ut ~ 10^^ Hz for BHBs and Vt ~ 10^^ Hz for AGNs, see e.g. 
WZOO. This corresponds to the value of Ut/uc from tens to 
hundreds, where fc = eB /2TimcC is the cyclotron frequency 
for the magnetic field of the strength B and rric and c are, 
respectively, the rest mass of the electron and the velocity of 
light. The ratio, Ut/uc, is larger in AGNs as a result of their 
magnetic field weaker than that in BHBs. The Lorentz fac- 
tor, 7t, of the electrons that emit most of their synchrotron 
radiation at Ut is of the order of 7t ~ {vt/vcY^'^ . The value 



Figure 2. Synchrotron spectra from a spherical source with 
the energy content of the non-thermal tail of (a) 5 = 0.19, (b) 
(5 = 10~*, and other plasma parameters as in our model of emis- 
sion from Cyg X-1 in Section ^.3| . The heavy short-dashed curve 
shows the optically-thin emission of the non-thermal electrons 
alone, (4/3)7til^ (47Tjj'), where j'j' is their emission coefficient. 
The heavy long-dashed curve shows the self-absorbed emission of 
the non-thermal electrons, iT(?B?S^^, where Sf} is their source 
function. The thin short-dashed and long-dashed cuves show the 
thermal (5 = 0) optically-thin and optically-thick emission, re- 
spectively. The heavy and thin solid curves give the corresponding 
emergent spectra, which change from optically thick to optically 
thin around the turnover frequency. 

of 7t is of the order of a few, again larger in the case of 
AGNs than of BHBs. 

First, let us consider a homogenous cloud of thermal 
plasma. The CS spectrum below the thermal turnover fre- 
quency, 1/*'^ can be approximated as a Rayleigh- Jeans spec- 
trum. Above , the synchrotron emission drops rapidly 
and the spectrum is often dominated by a power-law compo- 
nent resulting from Comptonization, which is cut off above 
an energy of ~ kT, where T is the electron temperature and 
k is the Boltzmann constant. 
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Now, let us assume that the electron distribution is 
given by equation (jl|). If 7nth < 7t, the number of electrons 
with 7 ~ 7t will increase (compared to the thermal case) 
as well as the absorption coefficient and in consequence the 
new turnover frequency, will be larger. Also the shape 
of the synchrotron spectrum is modified in this case. At low 
frequencies both emission and absorption are dominated by 
thermal electrons so that the optically thick spectrum is 
still the Rayleigh- Jeans one. However, at higher frequen- 
cies the emission of non-thermal electrons starts to domi- 
nate, while absorption is still dominated by thermal elec- 
trons [which is a consequence of the thermal distribution 
falling off much more steeply than the non-thermal one, see 
equation (^], and the optically-thick spectrum quickly in- 
creases. This leads to the regime in which both the non- 
thermal emission and absorption dominate, in which case 
the spectrum is proportional to the source function of non- 
thermal electrons. Finally, the spectrum becomes optically 



thin (c 



-{p-l)/2 



), turning over above 



This behaviour is shown in Fig. |^ where the resulting 
synchrotron spectra are calculated for two different energy 
contents of the non-thermal tail and, as a reference, for the 
purely thermal case. Note that the non-thermal tail has to be 
sufficiently strong for the optically thick spectrum to reach 
the non-thermal source function before it turns over. 

Both optically thick and optically thin parts of the syn- 
chrotron spectrum will undergo Comptonization. This pro- 
cess will thus be modified with respect to the thermal case 
due to two effects. Firstly, Compton scattering itself will in- 
volve both thermal electrons as well as non-thermal ones. 
The effect of non-thermal electrons on Comptonization will 
be most prominent at energies hv 3> kT (where h is the 
Planck constant), where a high-energy tail (above the ther- 
mal cut-off) will develop. Secondly, the seed-photon flux 
will be higher, so luminosity in the thermal Comptoniza- 
tion spectrum will also increase. Since 7t in optically-thin 
plasmas around black holes is expected to be well above the 
mean energy of the thermal electrons, even a non-thermal 
tail carrying a tiny fraction of the total electron energy, and 
thus only weakly modifying the shape of the Comptoniza- 
tion spectrum, can increase the turnover frequency and thus 
the luminosity produced by the CS process by a large factor. 

Then a question arises whether the presence of strong 
synchrotron self-absorption, acting as a thermalizing mech- 
anism for electrons of 7 7t (i.e. in the regime of optically 
thick emission), will allow for an electron distribution devi- 
ating substantially from a pure Maxwellian for 7 7t. This 
issue was addressed by Ghisellini et al. (1988), who have 
shown that an initial power-law distribution of injected elec- 
trons will develop (after the thermalization time scale equal 
to a few synchrotron cooling timescales, tl^'^, defined for the 
mean energy of the electrons radiating in the optically thick 
range of the spectrum) into a Maxwellian stationary distri- 
bution for 7 < 7t, provided the synchrotron cooling dom- 
inates and tl^'^ ^ tcsc (where tcsc is the escape time from 
the source). Assuming equipartition of magnetic field pres- 
sure with gas pressure, we find the latter criterion is always 
fulfilled in accretion fiows onto black holes, independent of 
the form of the injected distribution. 

However, when the inverse Compton process is dom- 
inant as a cooling mechanism, both 7nth and kT dimin- 
ish and a power-law tail starting below "ft is formed, see 



Ghisellini & Svensson (1990). A sufficient condition for the 
dominance of Compton cooling is a < 1, which is typically 
fulfilled in astrophysical accreting black holes. Even when 
Q > 1, Comptonization of external soft photons (e.g. from a 
cold accretion disk) can dominate electron cooling (see, e.g. 
Ghisellini et al. 1998). Thus, 7nth < 7t is a likely condition 
in those sources. 



4 SYNCHROTRON TURNOVER FREQUENCY 

Let fji, (7) be the synchrotron emission coefficient of a single 
electron averaged both over the electron velocity direction 
and the direction of emission. Then, the absorption coeffi- 
cient for an isotropic distribution of electrons in a chaotic 
magnetic field is (e.g. Ghisellini & Svensson 1991) is 



-1 



7 7 



no (7) 



7(7^ 



1)5 



d7- (5) 



We assume a homogenous spherical source of the radius, 
R, and the Thomson optical depth, tt, where rx = nctJT-R 
and (Tt is the Thomson cross section. Then, we can numer- 
ically calculate the absorption coefficient with equation (^ 
for the distribution of equation (hi) and then solve equation 
(^) for the turnover frequency, ^. This can be compared 
with the thermal turnover frequency, which for plasma pa- 
rameters typical for accreting black holes can be approxi- 
mated as (WZOO) 



th ^ ^ in3/^0 95 0.05 0.91 

ft ~ 7 X 10 B Tt i^c 



Hz. 



(6) 



On the other hand, we can obtain an analytical esti- 
mate, uf^, of the turnover frequency for the hybrid electron 
distribution by considering a purely power-law electron dis- 
tribution normalized to match the Maxwellian at 7nth- For 
the synchrotron absorption coefficient for power-law elec- 
trons (Rybicki & Lightman 1979), 



which, with equation IMj, leads to 



nt;h)7nth. 



(7) 



34+P2 t+PTC^+PUc 



7nth 



p-2 p-1 



6 + 159 GlrecrT^e7g~ '^ 
4 + 56 (TT 



where 



Gi 



(8) 



(9) 



Vc is the classical electron radius and F is Euler's gamma 
function. 

Provided p J; 3 (so that radiation from electrons with 
7 > 7f is negligible) and emission at i/^^^ is dominated by 
non-thermal electrons (i.e. 7f > 7t > 7nth, which condition 
can be checked a posteriori by comparing the source func- 
tion of power-law electrons to the Rayleigh- Jeans intensity) , 
t-""^ ~ uf^. On the other hand, uf^ < i^l^ when emission 
from thermal electrons is significant at u^^^ and, in general. 



j.f'^-maxfj.r,;^?') 



(10) 
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The accuracy of this approximation is typically ^ 10 per 
cent (provided 7t < -yt and p ^ 3). See WZOO for a discussion 
of approximations of i^l^. 

We can obtain the dependences of i/^^^/i^l^ on the 
plasma parameters by considering the ratio i/J^'/fJ'^ [given 
by approximations (^) and (^)]. First, at constant 5, this 
ratio decreases with increasing O since then 7nth increases 
with G [see equation (^)]. The dependence on tt is stronger 
for vf^ than for vl'^ as a result of much steeper optically 
thin synchrotron spectrum in the latter case. This leads to 
uf^/vl^ oc j ^jjg larger the optical depth, the 

larger the ratio jv^^ . In a hot accretion flow, we generally 
expect Tt increasing and B decreasing with increasing accre- 
tion rate, so the effect of a non-thermal tail (of a given form) 
on the value of the turnover frequency can be expected to 
be the largest for the most luminous sources. 



Since vflv't 



(2+p)/(4+p)-0.91 



the non-thermal in- 



crease of the turnover frequency grows with decreasing v^- 
This effect is due to the increase of 7t with decreasing 
[as a result of increasing ul^/uc, see equation (^], which 
leads to a much larger ratio of the number of non-thermal 
to thermal electrons at 7t for a given 7nth. Thus, the effect of 
a non-thermal tail with a given form will be much stronger 
in AGNs than in BHBs (due to much weaker magnetic field 
in the former objects). 

Fig. 1^ shows the values of Ut^^^/ul^ as a function of 
(7nth — l)/6 for three sets of thermal plasma parameters. 
In the first and second set, we assume parameters relevant 
for luminous BHBs, tt=1, 7? = 3 x 10^ cm, B = 10^ G, 
and = 0.1 and 0.2, respectively. In the third set, we as- 
sume parameters relevant to AGNs, tt~1, R = 3 x 10^^ 
cm, B = 10^ G and = 0.1. For each set, we consider 
two different slopes of the non-thermal tail, p = 3 and 4 
and 7f = 30. Additionally, we consider the cases of 7f — 10 
and 100 for p = 3 (for p > 3 the influence of the cut-off is 
much weaker). We see that even a weak non-thermal com- 
ponent, with (7nth — l)/0 ~ 12 (which corresponds to only 
~ 1 per cent of the total energy density of the electrons in 
the non-thermal tail) can lead to an increase of the turnover 
frequency by a factor of ~ 1.5-2 for BHBs and ~ 10 for 
AGNs. 

The dependence of the relative increase of the turnover 
frequency on p and 7f is relatively weak. This reflects the fact 
that the turnover frequency depends mostly on the number 
of electrons at 7t and since for small values of S, 7t is only 
slightly larger than 7nth, '^c(7t) depends weakly on either p 
or 7t. This suggests that the degree of the dependence on 
p and 7f should increase with increasing v^"^^ jv^ ^ in agree- 
ment with our results. It is important to note that weak 
dependence of Vt^^/ul^ on the shape of the non-thermal tail 
makes our results weakly dependent on details of the accel- 
eration mechanism. 
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5 COMPTONIZATION BY A HYBRID 
ELECTRON DISTRIBUTION 

5.1 The general case 

For Comptonization by a thermal plasma, we employ here 
a simple treatment of this process of Zdziarski (1985, 1986). 
In that approximation, the Comptonization spectrum above 



Figure 3. The increase of ft as a function of the electron energy 
where the non-thermal tail starts for various plasma parameters 
and forms of the tail. The upper axis shows the values of S for 
p = 3 and 7f = 30. Note that the relation between (7nth — 1)/© 
and (5 is slightly different for the curves with other values of p and 
7f , see Fig. ^. 
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Figure 4. The relation between tlie fraction of energy carried 
by non-thermal electrons and the ratio of Lj^,'/L^ for different 
plasma parameters and a = 0.7. 



the frequency, Vinj, at which soft seed photons are injected 
(which, in the case of the CS process, x vl ), can be ap- 
proximated as a sum of an e-folded power-law with energy 
index a and a Wien spectrum. 



dx 



r{2a + 3) 



-x/e 



(11) 



where x = hv/nriaC (hereafter all indices of x have the 
same meaning as those of v) and P^c is the volume- averaged 
scattering probability. Its integration yields the thermal- 
Compton luminosity, which, in the case of a spherical source, 
is 



r 1 



e 



+ 



6r(a)n, 



r(2a -I- 3) 



(12) 



where C is a constant depending on the flux of the injected 
seed photons (we calculate C in Section |5.2| for the case of 
synchrotron seed photons). 

The presence of non-thermal electrons modifies the 
Comptonization spectrum and since for the electron distri- 
butions we consider (i.e. rx ~ 1 and 5 <^ 1) the Thomson 
optical depth for scattering off non-thermal electrons is ^ 1, 
the resulting spectrum can be approximated as a convolution 
of the thermal Comptonization spectrum with a spectrum 
resulting from single scattering of photons off non-thermal 
electrons. Therefore, for a; 0, the Comptonization spec- 
trum becomes harder, while for x ^ Q>, above the ther- 
mal cut-off, a power-law tail in the spectrum develops. As a 
result, the overall luminosity produced by Comptonization 
increases. 

In calculations of the luminosity from Comptonization 
on non-thermal electrons, L^', where the Klein-Nishina ef- 
fect has to be taken into account, we use the approximation 
for the rate of energy change, A'y/At, of a single electron 
via Compton interaction with isotropic photons of energy 
density, [/ph, and mean energy, (a;), 



d7 _ 4 o-T (7^ - 1) t^ph 
At 3 meC 



1 



63 7(2^ 
To {x) 



(13) 



where the Klein-Nishina cross-section was approximated us- 
ing the first-order correction to the Thomson-limit cross- 
section (Rybicki & Lightman 1979). We assume that the 
photons undergoing scattering off non-thermal electrons are 
those from the thermal Comptonization spectrum above 
linj. Then 



3{AL'i^/Ax) (3/4 + rT/5) 



(14) 



where the factor (3/4 + rT/5) accounts for the change of 
the escape time due to scatterings in the source and is a 
matching formula between the optically thin case (where the 
escape time is 3i?/4c) and the optically thick one {ttR/5c, 
Sunyaev & Titarchuk 1980). We can further simplify the cal- 
culations assuming that the spectrum of photons undergoing 
non-thermal Comptonization is a pure power-law and then 
neglect scatterings in the Klein-Nishina limit (see below). 

Now, for each 7 we calculate Uph, {x) and {x^) as an 
integral over the power-law spectrum from a^inj up to the 
energy, Xmax, of the limit of the Klein-Nishina regime, which 
we obtain from the condition of 1 — 63/(107a;niax) = 0, in 
which the numerical coefficient is chosen for consistency with 
equation ([l^). 

Then the formula (|l^ is integrated over nc(7), from 
7nth to infinity, which leads to 



LtiC) 



167t 

~3~ 



7nth 



6 -use 

4 + 56 

a + l 
Tilth 



(!-t) 



7nth 



where equation (^ has been employed 
The total luminosity, Lg' = 



rpl 



(l-a)(2-a) 
(15) 



as well as the 



ratio / can now be directly obtained from expressions 



(^ and (|l5|). In the range of 0.4 ^ a 0.9, both the high- 
energy part of the power law component dominates the lumi- 
nosity as well as the Wien component can still be neglected. 
Then, we obtain (see WZOO) 



7-pi 
rth 



~ 1 

~ 3 

7nth 



/lOy-" 6 + 159 /3 Ti-\ 
163/ 4 + 59 U 5 / 

1 



2 - a 

7nth 



(16) 



We find the accuracy of this approximation is typically ^ 
30 per cent. Figure ^ shows that while the ratio Lq / Vq 
depends strongly on 5, the dependence on the non-thermal 
tail shape is very weak. 

In Fig. ^, we show as an example Comptonization spec- 
tra from a homogenous plasma cloud, produced with the 
Monte Carlo code of Gierlihski (2000) . We see that eis long 
as 5 < 0.1 the Comptonization spectrum at 3; 9 is hardly 
modified by the presence of non-thermal electrons and the 
onset of the power-law tail is no higher than an order of 
magnitude below the peak in the xL{x) spectrum. 

From the observational point of view, an important fea- 
ture measuring the deviation from purely thermal Comp- 
tonization spectrum is the fraction of the luminosity radi- 
ated in the high-energy tail beyond the thermal spectrum. 
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Figure 5. Monte Carlo Comptonization spectra of 10-eV black- 
body photons in a plasma of © = 0.1, rx = 2.4. The solid line 
represents the thermal spectrum, while the other lines correspond 
to the cases of power-law tails in the electron distribution of p = 3, 
7f = 30 and different values of S. 



The ratio L^/L^ is larger than this fraction as it incorpo- 
rates also the non-thermal luminosity at energies below the 
tail. This effect is stronger for smaller 5, when a larger frac- 
tion of Z/Q is hidden below the thermal peak. For example, in 
the spectra in Fig. ^ corresponding to S = 0.1 and 5 = 0.01 
the high-energy tail develops at ~ 300 keV and ~ 500 keV, 
respectively. Then L^'/Lj? is larger than he fraction of the 
luminosity in the tail by factors 3 and 5, respectively. Thus, 
J* provides only an upper limit to the fraction of the 



luminosity radiated in the high-energy tail. 



5.2 Comptonization of synchrotron photons 

As we have seen above, the presence of a non-thermal com- 
ponent in the thermal distribution of electrons in a plasma 
where synchrotron radiation is produced and then Comp- 
tonized leads to an increase of the luminosity available from 
this process in two ways. First, it increases the number of 
soft photons (both from optically thick emission up to x^^^ 
and from optically thin emission above a;"**") and second, it 
increases the efficiency of Comptonization, the former effect 
being significantly stronger. 

In the case of purely thermal synchrotron Comptoniza- 
tion, the normalization of the Comptonization spectrum can 
be assumed (Zdziarski 1985, 1986) to be proportional to the 
Rayleigh- Jeans flux at the turnover frequency, so that the 
constant C in the expression (p^), which we now denote as 
C"", reads. 



where 



(29)2 



1 + 10(29)2 



(17) 



(18) 



is a heuristic formula obtained by Zdziarski (1985) to match 
Monte Carlo simulation, and Ac is the Compton wavelength. 
The total luminosity from the CS process is then 



+ Lc(C"^) 



(19) 



where Lg^ is the luminosity in the thermal synchrotron spec- 
trum. 

If thermal and hybrid synchrotron photons experience 
the same average number of scatterings off thermal elec- 
trons before reaching a; ~ 9 (i.e. when Xt^^^/xf" < 9/2;^'' 
and the synchrotron spectrum is sufficiently narrow, which 
corresponds to p > 3 and which we assume hereafter in this 
section) , the ratio of the luminosity from Comptonization off 
thermal electrons in the hybrid case to that in the thermal 
case will equal the corresponding ratio of the synchrotron 
luminosities. Then, 



^nth 



r th 



(20) 



where L'g'^ is the hybrid synchrotron luminosity and C""' 
is the normalization constant for thermal Comptonization 
spectrum of hybrid synchrotron photons. 

To calculate the ratio ( |20[ ) we assume the optically thick 
hybrid synchrotron emission below a;""^ to be oc S'^^{x), the 
source function of power-law electrons, 

„3 1 , 



S^\x) 



1 G2 mcC^ 



231/2 Gi 



X3 



where 



G2 



r(^)r(^)r(2£^) 



f 3p+i9 ^ 
12 

r(^) 



The emission above is oc i 

3^(p + 4)G2 (<V^c)'^' 



(21) 

: 1. (22) 
(p-i)/2 g^jjjj then we have 



nth 



C 



14(p-3)Gi9(4Va:c 



(23) 



The above expressions were derived using the emission and 
absorption coefficients for power-law electrons (e.g. Rybicki 
& Lightman 1979). Since we neglected the cut-off in the 
electron distribution, the above relation holds for x^^^ /xc ^ 
jf, i.e. when the contribution from electrons with 7 > 75 to 
the emission at x^^^ can be neglected. 

The total luminosity of the CS process in the hybrid 
case is then 



Lnth 
CS 



Lnth I )rth/^nth\ 
S + ivc ) 



pi /^nth\ 



Lg'(C 



(24) 



When a < 1, we can neglect and is"^ equations ( p^ 
and (^) respectively, since most of the luminosity is then 
produced at much higher energies, a; 9. Then, the ratio 
Lq^/LJ^s can be calculated as 



-^cs 

rth 



^nth 

Tnir 



1 + 



(25) 



We can neglect the last term in the brackets if S is sufficiently 



small, in which case 



LS«^ 35(p + 4)(a;fV^c) 



7/2 



rth 



14(p-3)9(sJV^c)' 



C"M. We then have 



(26) 



Note that in calculations of C"**" we used the source 
function for power-law electrons. Thus, this approximation 
does not have the correct thermal limit when 7nth 2> 7t, i.e. 
when a;"*'' ft; a;*'' (see Fig. where a synchrotron spectrum 
for such a case is shown) . Note also that for harder power- law 
electron distributions {p ^ 3), for which the approximation 



© 2000 RAS, MNRAS 000, [[jl^ 



8 G. Wardzinski and A. A. Zdziarski 



10 




Figure 6. The amplification of tlie CS emission due to the pres- 
ence of nonthermal electrons, / L^^g{5), for models (see text) 
of NGC 4151, GX 339-4 and Cyg X-1 (dashed, solid and dotted 
curves, respectively). The horizontal lines correspond to the am- 
plification required to account for the observed X-ray luminosities 
and the corresponding values of 7nth are 2.61, 2.46, and 4.10 in 
the above 3 objects, respectively. The vertical line corresponds to 
S obtained from observations of the non-thermal tail in the spec- 
trum of Cyg X-1 by McConnell et al. (2000a), which corresponds 
to the X-ray luminosity exceeding that observed by a factor of 
~ 500. 

( po|) does not hold, the amplification of CS luminosity will 
be even larger than that predicted above. 

5.3 Applications to NGC 4151, GX 339-4 and 
Cyg X-1 

We now calculate the ratio LJ^g^/L^g as a function of the 
energy content of the electron power-law tail in three dif- 
ferent cases. They correspond to the models of thermal CS 
emission from a hot accretion disc with magnetic field in 
equipartition with the gas internal energy (dominated by 
hot ions) of WZOO for the Seyfert galaxy NGC 4151 and 
two BHBs in their hard states, GX 339-4 and Cyg X-1. The 
CS luminosities predicted by the thermal model (fig. 8a in 
WZOO) were orders of magnitude below the ones observed 
from NGC 4151 and GX 339-4, and less by a factor of ~ 3 
for Cyg X-1. 

For NGC 4151, we assume (see WZOO for details and 
references) R=1.2x lO" cm, B = 3.8 x 10^ G, a = 0.85, 
e = 0.1, p = 4, 7f = 100. For GX 339-4, the same values of 
G, p and 7f are adopted while R = 8.8 x lO'' cm, B = 1.5 x 
lO'' G, a = 0.75. In the case of Cyg X-1, R = 3x10"^ cm, B = 
6.3 x 10® G, a = 0.6, and we use the electron distribution 
obtained in McConnell et al. (2000a) by fitting the spectra 
from the COMPTEL and OSSE detectors aboard CGRO, 
i.e. O = 0.17, 7nth = 2.12, p = 4.5, at assumed 7f = 10'', 
which corresponds to a rather large value of 5 — 0.19. 



The results are shown in Fig. |6|. In the case of NGC 4151 
and GX 339-4, a non-thermal tail in the electron distribu- 
tion with S 10"'^ amplifies the power produced by the CS 
process sufficiently to reproduce the luminosities of those 
sources. Therefore, unlike the results obtained by WZOO 
from investigating the purely thermal CS process, in the case 
of electron distributions with relatively weak non-thermal 
tails the CS process appears to be capable of producing the 
observed X-ray spectra of accreting black holes, independent 
of the Eddington ratio or the black hole mass. 

On the other hand, the estimated thermal CS luminos- 
ity in the case of Cyg X-1 in the hot disc model is compa- 
rable to that observed from the source. Then, we find that 
the non-thermal tail observed by McConnell et al. (2000a) 
would produce the CS luminosity ~ 500 times higher than 
that actually observed (see Fig. Ba, where thermal and hy- 
brid synchrotron spectra are plotted for plasma parameters 
assumed in this model). The requirement that L^g cannot 
exceed the source luminosity constrains the magnetic field 
strength in the model to be B ^ 2.3 x 10^ G, which is 
~ 30 times smaller than the equipartition value estimated 
by WZOO. 

A similar constraint for Cyg X-1 can be obtained in a 
model of active coronal regions above the disc with magnetic 
field dissipation (see WZOO for details of the corresponding 
model of CS emission). Assuming A^ = 10 active regions of 
radius Rh = 3 x 10® cm and height Hh = 3 x 10^ cm, we 
obtain B ^ 3.4 x 10® G, much less than the field strength 
predicted assuming dissipation of magnetic field at a fraction 
of the Alfven speed, Bdiss = 5 x lO'^ G. Though varying N, 
Rh , Hh would change the upper limit on the field strength as 
well as the value of Bdiss, the relation between these quan- 
tities would remain similar. This suggests that the mecha- 
nism of dissipation of disc field in small active regions above 
the disc cannot be responsible for the energy release in the 
corona. This stems from the fact that -Baiss is the minimum 
field strength necessary for the magnetic field to dissipate 
enough energy so as to produce the observed luminosity. 
Therefore, either the model of active regions above a cold 
disc is ruled out for the hard state of Cyg X-1 or another 
mechanism of energy release must operate. We point out 
that if the CS process were to be negligible (as suggested by 
observations, see Section]^, the magnetic field would have 
to be even weaker, which would strengthen our conclusion 
even more. On the other hand, a weak magnetic field does 
not present analogous difficulties in the hot accretion disc 
model as the disc heating is directly by gravity. 

We then check how sensitive the above results are to 
details of the observed spectra of Cyg X-1. First, the spec- 
tral index of Cyg X-1 varies, and the value a = 0.6 assumed 
above is the hardest one and spectra with a — 0.7 have 
been also observed. Second, the relative normalization be- 
tween the COMPTEL and OSSE spectra remains relatively 
uncertain, and it may be possible that the former is less by 
~ 2 (McConnell et al. 2000b) than the one used in the calcu- 
lations above. This would then diminish the energy content 
in the electron non-thermal tail twice, which corresponds 
to the new value of 7nth ~ 2.31. Thus, we have repeated 
our calculations for a = 0.7 and 7nth = 2.31. The resulting 
constraints on the magnetic field strength in the hot disc 
and active regions models are then B ^ 5.8 X lO'^ G and 
B 6.8 X 10® G, respectively. These values are only about 
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twice the previous upper limits, which does not affect our 
conclusions above. 

We note that another constraint on the CS emission 
may be provided by extrapolating the X-ray power law spec- 
trum to lower energies (but ^ hvt), e.g. to the UV or optical 
ranges. That emission should not, in any case, exceed that 
observed. In the case of Cyg X-1, we have extrapolated the 
observed X-ray power-law spectrum of Cyg X-1 (Gierlihski 
et al. 1997) to v ~ 10^^ Hz and have found it is more than 
2 orders of magnitude below the UV spectrum of the com- 
panion supergiant as observed by lUE (Treves et al. 1980). 
Thus, those data are compatible with the CS origin of the 
X-ray spectrum in Cyg X-1. 

On the other hand, GX 339-4 is a low-mass X-ray bi- 
nary with a low upper limit on the emission of the com- 
panion star (e.g. Zdziarski et al. 1998). Then, simultaneous 
(as the object is highly variable, e.g. Corbet et al. 1987) 
observations in X-rays and at longer wavelengths can put 
constraints on the CS model in this case (e.g. Fabian et al. 
1982). 



6 CONCLUSIONS AND DISCUSSION 

We have studied the influence of non-thermal tails in the 
electron distribution on the Comptonized synchrotron emis- 
sion from plasmas of parameters typical for accreting black 
holes. In those plasmas, most of the self-absorbed thermal 
synchrotron emission is produced by electrons far in the tail 
of the Maxwellian distribution. As a result, even a weak non- 
thermal component far beyond the thermal peak can greatly 
enhance the synchrotron emission. This emission provides 
seed photons for Comptonization, which is also correspond- 
ingly enhanced. However, the shape the Comptonization 
spectrum remains only weakly modifled, with the main sig- 
nature being a relatively weak high-energy tail. 

We have shown that this effect becomes stronger with 
the decreasing plasma temperature and increasing optical 
depth, which typically correspond to an increasing Edding- 
ton ratio. It is also much stronger in the case of AGNs than 
of accreting stellar black holes, as a result of weaker mag- 
netic field in the former objects. The CS luminosity can be 
amplified by a factor ~ 10^ and ~ 10'^, respectively, in the 
case of only ~ 1 per cent of the electron energy in the non- 
thermal tail. 

Then, the Comptonized synchrotron omission becomes 
energetically capable of producing the observed luminosities 
of accreting black holes independently of their Eddington 
ratio or the black-hole mass. This conclusion is qualitatively 
different from that obtained in the case of purely thermal 
plasma. 

It is then possible to obtain upper limits on the mag- 
netic field strength in optically thin plasmas around black 
holes by constraining the parameters of the non-thermal 
component in the electron distribution from observations 
of the high-energy tail in Comptonization spectra. Based 
on the CGRO data for Cyg X-1 in the hard state, we have 
obtained the upper limits on the field strength at least an 
order of magnitude below both the value corresponding to 
equipartition (in the model of a two-temperature accretion 
disc) and the minimum value required by dissipation of mag- 
netic fields (in the model of active coronal regions). There- 



fore, the latter model appears to be ruled out for the hard 
state of Cyg X-1. 

We note that the presence of a strong correlation be- 
tween reflection strength and the X-ray spectral index, ob- 
served in Cyg X-1 (Gilfanov et al. 1999), suggests that the 
CS emission is a negligible source of photons for Comp- 
tonization and thus the magnetic field strength must remain 
significantly below the upper limits derived above. The same 
correlation seen in other objects (Zdziarski et al. 1999; Gil- 
fanov et al. 2000) implies that either the high-energy tails 
or the magnetic fields are weak enough for the CS process 
not to dominate the energy output. 

In this work, we have assumed, for the sake of simplic- 
ity and compatibility with studies of, e.g., McConnell et al. 
(2000a), the electron distribution of the non-thermal elec- 
tron tails to be cx 7"'' at any value of 7. However, when 
the power-law tail extends down to non-relativistic energies, 
7nth ~ 1, the distributions expected from acceleration pro- 
cesses are power laws in either the kinetic energy or the 
momentum, e.g. oc (7 — 1)"*" for the former. Clearly, the 
smaller 7nth (i.e. the larger 5), the larger the differences 
between those distributions. Then, the actual shift of the 
turnover frequency at a given value of 7nth would be some- 
how smaller. Also, the relation between 7nth and 5 would 
change. Nevertheless, the above effects would modify our re- 
sults only quantitatively without affecting our conclusions. 

The influence of a weak, non-thermal, component in the 
electron distribution on radiation spectra has also been in- 
dependently considered by Ozel, Psaltis & Narayan (2000). 
Their study has been devoted to the case of advection- 
dominated accretion flows (ADAF), and their presented 
spectra are integrated over all radii of the flow. In contrast 
to our results, they do not note any shift in the turnover 
frequency at the peak of the synchrotron spectrum. This 
appears to result from their study being constrained to the 
ADAF model at low accretion rates, in which case electrons 
in the innermost parts of the accretion flow reach rather 
high temperatures (at which even thermal electrons reach 
relatively high Lorentz factors). Those temperatures are sig- 
niflcantly higher than those we consider based on obser- 
vational data from luminous accreting black holes. On the 
other hand, Ozel et al. (2000) find a signiflcant excess non- 
thermal emission at frequencies well below the peak of the 
integrated synchrotron spectrum. This excess results from 
radiation emitted at large radii, where the electron tem- 
perature is low, and where the synchrotron emission (both 
optically thick and optically thin) around the local turnover 
frequency can be strongly amplified by the presence of an 
electron tail, as illustrated in Fig. above. 
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